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INTRODUCTION 

T*V\o  Ur»r*l/  nraeonfurl  Wnfo  {  e  t*n  1  n  f  ►/>  liar  J  mia  nr*rtK  1  amA  onnnarna/l  i.»  V» 

the  propagation  of  shock  waves,  mainly  Imploding  spherical  and  cylindrical 
shock  waves,  in  water  and  (in  part)  in  gases. 

In  Section  1  the  relevant  equation  of  state  for  water  is  presented  and 
the  following  section  deals  with  the  consequences  of  this  equation  on  shock 
wave  transitions. 

Section  2  is  concerned  with  the  self-similar  solution  valid  in  the  final 
stages  of  collapse,  where  the  shock  may  be  considered  to  be  infinitely  strong 
so  that  the  pressure  ahead,  and  more  Important,  a  pressure  term  in  the  equa¬ 
tion  of  state  mry  be  neglected.  Numerical  results  are  given  for  the  expo¬ 
nent  in  the  power  law  for  the  shock  speed.  In  the  following  section  a 
perturbation  on  this  solution,  linearized  in  terms  of  the  previously  neg¬ 
lected  pressure  terms,  is  evaluated.  In  effect  this  extends  the  range  of 
validity  of  the  self-similar  solution  further  away  from  the  collapse  point 
or  axis.  The  work  of  these  two  sections  is  closely  related  to  the  propaga¬ 
tion  of  shocks  and  detonations  in  gases  (Welsh,  1966),  where  a  much  fuller 
account  is  presented. 

An  approximate  solution  for  the  shock  motion  from  the  time  of  its 
initiation  to  its  final  collapse  is  given  in  Section  5,  the  method  employed 
being  similar  to  that  used  by  previous  authors  (Chester,  et.al.)  for  the 
motion  of  a  shock  wave  in  a  non-uniform  tube. 

The  final  section  is  devoted  to  the  exact,  linearized  solution  for  the 
Initial  motion  due  to  the  release  of  a  spherical  diaphragm  separating  two 
uniform  regions  at  different  pressures.  The  analysis  is  applicable  to 
i)  spherical  or  cylindrical  geometry,  ii)  high  pressure  in  the  interior  or 
exterior,  and  ill)  any  combination  of  gas  and  water  for  the  two  media. 
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1 .  EQUATION  OF  STATE  FOR  WATER 

In  the  classical  theory  of  hydrodynamics  the  compressibility  of  water 
is  neglected,  due  to  the  relatively  low  order  of  magnitude  of  the  pressure 
variations  involved.  However,  when  water  is  subjected  to  extremely  high 
pressures,  as,  for  instance,  in  an  explosion,  the  compressibility  can  no 
longer  be  neglected.  In  such  a  situation  the  Tult  equation  of  state  is 
known  to  give  an  adequate  representation  of  the  compressibility  (Cole,  1948). 
This  equation  relating  pressure,  entropy,  and  density  can  be  expressed  as 


p  +  B(s)  -  B(s)  ^  ^ 


(1.1) 


B(s)  is  a  slowly  varying  function  of  entropy  s  and  has  the  dimensions  of 
pressure,  p'  is  the  (finite)  value  of  the  density  for  which  the  pressure 
p  is  zero.  7  is  dimensionless  and  approximately  constant.  With  B(s) 
taken  to  have  the  constant  value  3.047  kilobars  and  7  taken  ns  7.15  the 
above  equation  of  state  gives  a  good  approximation  to  the  behavior  of  water 
at  pressures  up  to  100  kilobars.  Here  we  shall  take  7  ■  7. 

Defining  the  sound  speed  of  water  c  by 


■cn 


we  obtain  from  (1.1) 


.  !Le±b^2L 

p 

For  a  perfect  gas  the  equation  of  state  is 

s/c  v 
'  v  7 
p  «  e  p' 


(1.2) 


(1.3) 


. 
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If  we  assume  that  Che  function  B(s)  in  (1.1)  has  the  constant  value 
B  and  that  p  is  sufficiently  large  that  wc  may  neglect  B  on  the  left 
hand  side  of  (1.1),  then  the  equation  of  state  for  water  is  simply 

P  *  P7 

which  is  the  equation  of  state  of  (1.3)  for  a  gas  when  the  entropy  Is 
uniform,  i.e.,  the  motion  is  homontropic.  Also  the  sound  speed  (1.2) 
reduces  to 

c2  -  12- 

P 

so  that  the  motion  of  water  at  sufficiently  high  pressures  may  be  treated 
as  that  of  a  gas  with  7*7,  neglecting  the  effect  of  any  variations  in 
entropy.  However,  p  must  not  be  so  large  that  (1.1)  is  not  valid. 

These  simplifications  will  be  used  in  evaluating  the  self-similar  solution 
for  the  final  collapse  of  spherical  shock  wave. 

TVi** i  s  are  evidently  three  distinct  ranges  of  pressure.  The  first  for 
p  «  B  corresponds  to  incompressible  flow.  The  present  work  will  be  con¬ 
cerned  with  the  remaining  two,  i.e.,  p  *>  B  and  p  »  B.  B  will  be 
assumed  constant. 

We  shall  later  require  an  expression  for  cue  specific  internal  energy 


of  water,  defined  by 


4 


.  f 

J 


n  d 


f  ^ 

\p  / 


— - - -f-  £ 

7(7-1)  p 


•  •>/ 


2.  SHOCK  TRANSIT IONS 

The  propagation  of  a  shock  wave  in  water  is  governed  by  the  mechanical 
laws  of  conservation  of  mass,  momentum,  and  energy.  The  shock  is  assumed 
to  be  a  discontinuity  in  the  flow  variables,  viscosity  and  heat  conduction 
being  neglected.  Also  the  relations  between  the  flow  variables  on  either 
side  of  a  plane  shock  will  be  used  for  curved  shocks,  under  the  assumption 
that  the  width  of  the  shock  is  negligible  in  comparison  with  its  radius  of 
curvature. 

With  the  subscript  o  denoting  the  variables  ahead  of  the  normal 
shock,  the  variables  behind  it  are  given  by  the  three  conservation  laws 


p(U-u)  -  P  (U-u) 
o  o 


p  +  p(U-u)2  •  pQ  +  PQ(U*U0)2 


1  ,2  ,  £  ,  _ c_  ,  B 

2  (U-u)  +  p  +  7(7-1)  +  p 


1  ...  Y2  ^  Po  ,  Co  .  B_ 

2  (U‘V  +  p  +  7(7-lT  +  p 

O  0 


2  <“-“>2  •  1  <u-"„>2  +  h 


Also  the  second  may  be  written 


or 


2  \2  U-u  ,  2  ,  .,.,2. 

c  +  7  (U-u)  -  jTjj-  (cq  +  7U  ) 

c 
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so  that  D  does  not  appear  if  the  equations  arc  expressed  in  terms  of  u, 

c,  p,  in  which  case  they  coincide  with  those  for  a  gas.  We  shall  take 

u  3  0 , 
o 

The  solution  in  terms  of  the  pressure  ratio  z  *  p/pQ  is 


£_ 

P_ 


&&  +  (7-1)  +  (7+1)2 

Pq _  _ 


27B 


+  (7+0  +  (7-l)z 


21  +  2sl  !°  +  m  !a 

P„  2  0O  2 
O  0  0 


2  p0  (z-n* 
u  -  —  •  „„  *  -  r 


p  21  .  Zrl  7+1 

o  i—  +  4-  — —  Z 

p  z  z 


(7+l)2c2  -  jjp  (72+47-l)  +  472  ^  +  7(7Z-1)  ^  2 

kO  O  o 


.  12_l2_iL2 

p„ 


<pc+b) 


7B  +  ^  po  +  ^  Poz 


The  special  case  B  ■  0  corresponds  to  a  gas. 

In  terms  of  the  Mach  number  M  -  u/c  the  variables  behind  the  shock 

o 

arc 


Q_  „  (7+l)M 

Po  (7-l)MZ+2 
2c 


u  * 


£("-0  . 


£_  .  27M..-7+1  ,  27.  B_  ,  v 

P„  7+1  7+1  p  l) 

A  1 


(7+1) 


2  c 


[27M2-(7-l)H(7-l)M2+2] 


M 
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The  shock  Mach  number  M  in  terms  of  the:  pressure  ratio  is 


Since  Pq/B  ■  1/300,  if  the  water  is  initially  at  atmospheric  pressure, 

3 

the  factor  multiplying  2-1  1j  approximately  3.10  .  Hence  M  does  not 
differ  greatly  from  unity  unless  z  iu  of  corresponding  order  of  magnitude. 
In  the  limit  as  z  1,  M  1  and 

/V(p  +  B) 

U  -t  c  ■  /  -  » 

O  N  p 


so  that  a  very  weak  shock,  which  is  approx imately  sonic,  travels  at  a  much 
greater  speed  than  a  weak  shock  in  a  gas,  but  it  requires  a  much  greater 
change  in  pressure  ratio  to  increase  the  spaed  of  the  shock  wave  in  water 
appreciably . 


Strong  Shock  Relations 

A  useful  simplification  of  the  conservation  equations  for  a  shock  wave 
in  a  gas  is  obtained  by  assuming  that  the  shock  is  'strong'  so  that  the 
pressure  ahead  may  be  neglected  in  comparison  with  the  pressure  behind. 

This  simplification  is  essential  in  obtaining  self-similar  solutions  in¬ 
volving  shock  waves  whose  speed  is  non-constant.  In  the  case  of  water,, 
Identical  strong  s.hock  relations  may  be  obtained,  but  under  the  more  severe 
restriction  that  not  only  is  p  »  but  also  p  »  B.  These  are 

Bl.  m  Z±i 

7-1 


u 


2 

7+1 


U 


2 

c 


27(7-1).  y2 
(7+1) 2 
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3 .  SELT-SIMILAR  CO  LAPSING  SHOCKS 

The  sc J f -similar  solution  for  the  final  stages  of  a  collapsing  spher¬ 
ical  shock  wave  in  a  gas  is  well  known  (Guderlcy,  19A2).  A  similarity 
hypothesis  is  obtained  by  using  the  strong  shock  relations,  so  tTlat  at 
any  instant  the  flow  in  determined  only  by  the  radius  of  the  shock  and 
the  initial  density  of  the  gas,  apart  from  scaling  factors. 

The  equations  governing  the  symmetric  motion  of  a  perfect,  inviscid, 
non-heat-conducting  gas  arc 

<)  ,  . .  .  .  3  !•>" 


(u±kc)  +  (u±c)  |g-  (uike)  -  +  ^  +  f-  gf 

5t  +  u  3r  “  0 


(3.1) 


where  k  *  2/(7-l),  ?  »  s/[cv(7-l) J ,  where  s  is  specific  entropy,  and 
j  ■  1  for  cylindrical  symmetry,  2  for  spherical  symmetry.  The  appropriate 
equations  for  water  are  the  above,  neglecting  entropy  variations,  and  so 
are 


It  (u±kc)  +  (u±c)  (u±kr.)  -  +  — — 


(3.2) 


No  dimensional  constants  appear  in  these  equations  so  that  ihe  shock 
speed  U  is  dependent  only  on  the  radius  coordinate  X  of  the  shock 
front,  the  time  t,  and  the  density  P  .  Hence, 

A. 


U  * 


i.e. , 


t  «  >Ma 


or 


u  . 


and,  by  suitable  choice  of  length  scale, 

.  xl-(l/0£) 


U 


-  a  \ 


i/a 


Due  to  the  form  of  the  strong  shock  relations,  the  flow  is  self- 
similar  with  similarity  variable 


i 


a  r1//q 


and  the  dimensionless  fluid  velocity  and  sound  speed  r,  s  defined  by 


r  s  —  u 

a  R 


8  “  a  R  c 


are  of  the  form 


r  *»  r(4) 


so  that  the  equations  of  motion  are 


i  4r  (r±ks)  =  - — 4- 

d5  1-r+s 


s(l) 


where 


b+  =  (l-a(r±s)}(r±ks)  +  jars 


We  can  eliminate  4  from  these  equations  and  consider  r  as  a  function 
of  s  given  by 


1  dl 

k  ds 


(l-r+s)b+  +  (l-r-s)b 
(l-r+s)b+  -  (l-r-s)b 


I  .  r(r-l)far-l)  +  sZ[kq-a)-a(i+l)rl 
(1-r) [k-a(j+k+l)r]  +  r(l-ar)-tos^ 


(3.3) 


The  parameter  a  is  as  yet  unknown  and  is  determined,  as  in  the 


Gudcrlcy  solution,  by  the  condition  that  the  flow  is  regular  on  the 


negative  characteristic  which  reaches  the  center  of  symmetry  at  the  same 
instant,  t  *  0,  as  the  shock  wave.  On  this  line  r  *  l+s  and  the 
regularity  condition  is  simply  b_  ■  0  or  r  •  1+8.  The  solution  has 
to  satisfy  this  regularity  condition,  at  s  -  sq  say,  arid  the  conserve 
tion  equations  at  the  shock,  where 


7+1 


If  we  seek  a  regular  expansion  for  r(s)  about  s  »  bq,  say 

r  »  r  +  r.(s-s  )  +....+  r  (s-s  )n  +  .  .  .  . 
o  i  o  no 

then  rQ,  r^  are  the  solutions  of  quadratic  equations.  The  roots  for  tq 

are  real  in  some  range  CtR  <  a  <  1  say,  and  the  roots  for  r^  will  also 

be  real  unless  this  singular  point  (ro,SQ)  *n  r-s  plane  is  a 

spiral.  Thus,  there  are  in  general  four  solutions  which  are  regular  on 

the  limiting  negative  characteristic.  The  selection  of  the  solution  is 

decided  partly  by  trial  and  error  by  computation  and  partly  by  the  behavior 

of  the  integral  curves  of  the  equation  (3.3)  to  be  discussed  later. 

In  order  to  integrate  the  equation  (3.3)  a  similar  procedure  to  that 

used  previously  (Welsh,  1966)  was  adopted.  In  this  the  regular  series 

solution  is  developed  iteratively,  without  evaluation  of  the  actual 

coefficients,  from  the  1.  n.  c.  s  -  s  to  the  shock  s  «  -E,  and  the 

o 

discrepancy  there  noted.  The  actual  solution  and  the  appropriate  value 
of  0:  may  then  be  evaluated  by  trial  and  error. 

Exact  results  were  found  by  this  method  and  are  given  in  Table  I  for 
twelve  cases,  j  =  1  and  j  =  2  and  7  =  6/5,  7/5,  5/3,  3,  5,  7. 

In  the  case  of  a  gas  the  extremely  simple,  approximate  evaluation 
of  a  given  by  Whitham  (1953)  gives  remarkable  accuracy.  The  method 


10 

Is  applicable  here  and  consists  of  applying  the  characteristic  condition 
to  he  satisfied  by  the  variables  on  the  1.  n.  c.  to  the  values  at  the 
shock.  The  character is  tic  condition  is 


d(u-kc)  ■  dt 


on 


dR 

dt 


u-c 


and  on  the  shock 


u 


_  _2_ 
7+1 


c  « 


E  r1'*1/0) 


and  substitution  of  these  into  the  characteristic  condition  gives 

l.k  m  _ -JPg  2 

a  (D+E)  (D+kE)  •  D  *  7+1  (3.4) 

which  is  tabulated  as  along  with  the  exact  results.  The  relation  be¬ 

tween  aw  and  the  exact  value  is  apparently  more  systematic  than  for  the 
case  of  a  gas,  in  which  entropy  variations  occur.  In  all  cases  lies 

between  aR  and  the  exact  values  which  are  themselves  never  far  apart, 

and  the  accuracy  increases  with  decreasing  7.  For  7*1,  a  *1  as 

w 

is  the  exact  value.  For  this  special  case  the  shock  is  in  fact  a  character¬ 
istic.  Also,  the  Internal  energy  is  infinite,  causing  the  collapse  to  be 
uniform. 

The  Integral  Curves  and  the  Reflected  Shock 

The  equations  (3.3)  apply  to  the  homentropic  self-similar  flow  of  a 
gas  or  the  self-  similar  flow  of  water,  under  the  present  assumptions.  They 
are  identical  to  those  used  by  Hunter  (1960,1963)  for  the  collapse  of  an 
empty  spherical  cavity  in  water,  apart  from  his  use  of  the  dimensionless 
square  of  the  sound  speed.  We  shall  now  consider  the  behavior  of  the 
integral  curves  of  (3.3). 
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Due  to  rho  similarity  of  the  flow  the  variables  are  functions  of  4 
only.  The  physical  R-i.  plane  may  be  considered  to  be  made  of  linea  (para¬ 
bolas)  4  “  constant.  The  shock  is  4  ■  -1  and  5  <  "1  is  the  stationary 
region  ahead  of  the  shock.  The  speed  of  a  line 

„  t 

i  -  —  i7ri  “  co1181®06 
a  R  ' 

is 


Hence  the  nondimensionalization  of  the  variables  u,c  to  r,s  is  with 
respect  to  the  speed  of  the  line  5  ■  constant,  and  such  a  path  in  the 
R-t  plane  gives  rise  to  a  point  in  the.  r-s  plane.  The  shock  front  trans¬ 
forms  into  the  point 


r  ■ 


2 

7+1 


8  «  ■Jjw-ji 

8  7+1 


Negative  values  of  s  correspond  to  negative  values  of  t,  i.e.,  to  the 
imploding  phase,  and  positive  values  of  b  describe  the  expanding  phase. 

The  positive  and  negative  characteristics  of  the  flow  each  correspond  to 
a  point  on  the  lines  r  “  1  +  s,  respectively,  on  which  dr/dS,  ds/d{ 
are  singular  in  general.  Any  physical  solution  must  hove  4  increasing 
throughout  the  associated  integral  curve  in  the  r-s  plane.  In  general 
the  direction  of  increasing  time  on  the  integral  curves  in  the  r-s  plane 
changes  on  crossing  r  -  1  ±  s.  In  other  words,  a  physical  solution  cannot 
in  general  cross  these  lines.  However,  the  appearance  of  singular  points 
of  the  equation  (3.3)  for  dr/ds  on  r  -  1  ±  s  provide  the  means  of 
crossing  these  lines.  In  effect,  the  singular  points  cause  another  change 
of  direction,  ';«ncelling  that  appearing  already.  The  regularity  condition 
mentioned  previously  is  in  fact  the  condition  for  a  regular  solution  through 
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one  of  these  singular  points.  For  a  <  a„  these  points  do  not  appear,  for 
a  -  0^  they  coincide,  and  for  a  >  they  are  real  and  distinct.  If 
is  real  then  there  is  the  possibility  of  a  regular  solution  through  the 
singular  point,  which  is  either  a  node  or  a  saddle.  If  r^  is  imaginary 
then  the  point  must  be  a  spiral. 

The  finite  singular  points  of  (3.3)  are  as  follows: 

P4<°-0>-  P5  (s*  ao+k+l)  )  ’  P6(°'a) 


and  three  other  points  which  are  the  mirror  images  of  P^,  P^*  P^  in  the 

r-axis.  Both  sq±  <  0,  and  S  is  the  negative  root  of 


s2  «  g.CL:g.  -O 

k  a 


r  “  a(j+k+i) 


P^  is  always  a  degenerate  note  and  P^  is  a  saddle  point  with  asymptotes 

in  the  direction  of  the  r,s  axes.  is  always  a  node  and  is  degenerate 

if  k  -  j  +1  (i.e.,  7  ■  5/3  spherical,  7  •  2  cylindrical).  For  k>  j+1 

the  major  axis,  i.e.,  the  direction  of  all  the  integral  curves  except  for 

one  at  right  angles,  is  in  the  r-direction  and  for  k  <  j+1  in  the 
s-direction.  For  P^  there  are  two  possibilities.  When  it  is  above 
r  »  1+s  it  is  a  node  and  a  saddle  when  below,  the  transition  occurring 
through  its  merging  with  or  **3  on  r  “  1+s.  Except  in  the  special 

case  k  -  j+1  the  points  P^P^  seParate  into  a  node  and  saddle  or  saddle 
and  node,  respectively,  for  a  slightly  greater  than  d^.  When  they 
coincide  with  OL  ■  they  form  a  confluence  of  a  node  and  saddle.  This 
type  of  singular  point  is  not  a  simple  one.  In  Investigating  the  behavior 
of  the  curves  near  such  a  point  it  is  necessary  to  include  second  order 
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terms  In  the  numerator  and  denominator  as  the  discriminant  of  the  coeffic¬ 
ients  cf  the  flint  uuier  Lerms  vanisnes,  and  to  tirst  order  the  solution 
obtained  is  a  straight  line.  For  such  a  second  order  (singular  point  the 
significant  variations  occur  on  parabolas.  For  a  first  order  node  or 
saddle  the  field  of  curves  may  be  determined  by  finding  the  asymptotes 
and  the  lines  on  which  the  curves  have  zero  and  infinite  slope.  For  a 
second  order  singular  point  the  latter  two  lines  and  the  asymptote  between 
them  arc  parabolas  and  the  relative  positions  of  the  curves  are  the  same 
on  cither  side  of  the  point  (in  the  first  order  case  the  positions  are  re¬ 
versed)  resulting  in  a  naddle,  This  situation  is  exhibited  graphically 
in  Ffg.  1.  In  the  special  case  k  ■  j+1  there  is  a  triple  coincidence 
oi  P2,  and  P5  when  a  «  a^,  resulting  in  a  third  order  singular 

point,  so  that  the  significant  variations  take  place  along  cubic  curves. 

The  position  of  these  differs  on  either  side  of  the  point  and  the  situation 
is  geometrically  similar  to  the  first  order  case,  except  for  the  scale  in¬ 
volved.  Such  a  point  may  he  either  a  third  order  node  or  saddle.  Evidently 
higher  order  singular  points  will  be  geometrically  similar  to  the  first 
order  type  of  singular  point  (node  or  saddle)  or  the  second  order  type 
(naddle),  depending  on  whether  the  order  is  odd  or  even. 

The  values  of  a  for  which  these  transitions  just  discussed  occur  are 
tabulated  in  Table  I.  For  k  <  j+1  the  point  P^,  which  is  a  node  when 
it  first  appears,  changes  to  a  spiral,  in  which  case  no  regular  solution 
through  it  is  possible.  However,  the  actual  solution  was  found  to  pass 
through  in  all  cases.  Through  a  first  order  node  or  saddle  there  are 

two  regular  solutions.  In  the  case  of  a  saddle  these  are  the  asymptotes, 
which  are  the  only  two  solutions  passing  through  the  actual  point.  For  a 
node  one  of  the  two  is  the  minor  axis,  i.e.,  the  one  solution  not  having 


t 
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Che  same  alope  as  all  the  others  at  Che  point,  but  the  other  la  one  of  the 
remaining  curvaa  *n<|  her.ee  is  loiia  tllscinguishabie.  In  all  cases  the  actual 
solution  mss  through  the  major  axis  of  P^.  and  this  was  found  to  be  the 
only  possible  solution.  No  transition  from  P2  to  occurs  as  it  does 

In  the  case  of  a  gas  (Welsh,  1966).  Also  the  argument  given  in  that  pub* 
lication  that  no  uniform  collapse  (Ct  -  1)  is  possible  applies  here  to  the 
case  of  water,  since  the  flow  behind  a  uniform  shock  Is  homentropic  and  the 
equations  of  motion  for  gas  and  water  are  Identical.  This  contrasts  with 
the  situation  found  by  Hunter  (1960,1963)  for  cavitation  in  water.  An 
accelerating  collapao  exists  in  a  certain  range  of  7,  a  uniform  one  in  a 
different  range,  the  latter  being  singular  at  the  cavity  front  for  certain 
discrete  values  of  7 . 

The  appropriate  integral  curve  representing  the  actual  solution  in  a 
given  physical  solution  starts  at  the  shock  point  (D,-E),  crosses 
r  -  l+s  at  P2,  and  then  runs  into  P^(0,0),  with  time  increasing  through¬ 
out  this  path.  corresponds  to  the  line  t  *  0  in  the  physical  plane, 

and  the  solution  must  continue  beyond  P^.  The  integral  curves  are  symmetric 
about  r  *  0,  except  for  a  change  in  direction  of  increasing  time  due  to 

the  change  in  sign  of  s.  Thus  the  Integral  curve  passes  through  P^ 

smoothly,  the  points  of  the  curve  beyond  P^  corresponding  to  lines  £  >  0, 
i.e.,  expanding  lines  in  the  flow  behind  the  shock.  However,  the  curve 
heads  towards  the  singular  line  r  ■  1-s,  corresponding  to  a  positive 
characteristic  in  the  flow,  which  indicates  the  presence  of  a  reflected 
shock  wave  traveling  along  £  ■  constant  >  0.  The  solution  is  represented 
by  a  portion  of  this  curve  in  the  region  s  >  0,  r  <  1-8,  but  Jumps  to 
some  other  curve  above  r  ■  1-s.  A  condition  on  the  flow  behind  this  re¬ 
flected  shock  is  that  the  fluid  velocity  be  zero  at  the  center,  i.e., 
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u  =  0  on  R  ■  0.  Hence  r  is  finite  here,  but  since  c  i*  ~ 

is  infinite.  It  will  now  be  seen  that  this  Information  is  sufficient  to 
determine  the  integral  curve  which  satisfies  these  conditions  at  R  •  0, 
by  examining  the  singular  points  of  (3.3)  at  infinity. 

Here  we  are  interested  in  the  region  s  >  0,  There  are  three  singular 
points  on  s  -  +«  at  r  «  ±»,  ^irst  two  are  nodes  (half- 

nodes)  and  the  one  between  them  is  a  saddle  (half-saddle),  so  that  all 
curves,  save  one,  approach  either  r  »  ±o»  as  s  -»  «.  The  asymptote  through 
the  saddle  point  is  the  only  curve  which  has  a  finite  value  of  r  for 
s  «  +®  and  so  represents  the  flow  behind  the  reflected  shock.  The  Jump 
from  this  curve  to  the  one  describing  the  flow  ahead  of  the  reflected  shock 
Is  determined  by  the  conservation  equations.  In  this  way  the  value  of  £ 
on  the  reflected  shock  can  be  evaluated.  This  has  not  been  done  here. 

This  type  of  reflection  occurs  when  the  shock  speed  is  singular  and 
its  curvature  is  infinite.  If  there  is  a  symmetrically  placed  solid  spher¬ 
ical  or  cylindrical  boundary  causing  the  shock  to  be  reflected  before  its 
final  collapse  then  this  type  of  reflection  is  plane  in  nature,  provided 
the  radius  of  the  boundary  is  large  in  comparison  with  the  shock  width. 

4.  THE  PERTURBATION  OF  THE  SELF-SIMILAR  SOLUTION 

The  self-similar  solution  is  valid  for  p  »  B,  which  is  fairly  re¬ 
strictive,  and  presumably  means  that  it  is  valid  for  a  collapsing  shock 
wave  in  water  for  a  much  smaller  neighborhood  of  the  collapse  point  than 
in  the  case  of  a  gas.  Here  we  shall  consider  the  perturbation  on  the  self- 
similar  solution  due  to  taking  B  into  account  to  first  order.  The  pertur¬ 
bations  are  analogous  to  those  due  to  the  counter-pressure  ahead  of  a 

» 

collapsing  shock  wave  in  a  gas  (Welsh,  1966),  the  perturbations  being  of 
relative  order  R"2+(2/a>, 
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The  form  of  the  solution  is  now 

u  «  Lttl  Rl-(l/a)  +  rill  R-l+(l/a) 
c  .  Sill  Rl-(l/«)  +  spi  p-l+(l/a) 


and  the  equations  governing  these  perturbations  are 


£(l-r+s)  —  (r±ks) 

5(l-r+s)  (r±k?) 

where 

a±  -  +  ja(ra+r?)  +  (l-r*a)  (rtkiF)  +  (1-a)  (k-1) (rs-rs) 

+  -*-*■ b  -  (r±s)(r±ks)  , 

1-t+s 

and  £  may  be  eliminated  as  previously  by  considering  r,  r,  s  to  be 
functions  of  a 


(l-r-s)b(r ’+k)  «  (l-r4-s)b+(r'-k) 

a  (r'+k)  ■  b,(r'+k) 

t  *r 

a(r'-k)  -  b_(r,-k) 


where  '  =  d/4s 


\ 

- 

J 


(4.1) 


As  is  the  case  in  the  basic  self-similar  solution  the  boundary  condi¬ 
tions  are  obtained  from  the  regularity  condition  at  s  ■  sq  and  the  con¬ 
servation  equations  at  s  ■  -E.  Since  the  displacement  of  these  lines  from 
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their  positions  in  the  self-similar  flow  gives  rise  to  terms  of  the  same 
order  as  the  perturbations,  it  is  necessary  to  take  into  account  thee*-  dis¬ 
placements  wticn  applying  boundary  conditions  on  these  lines. 

In  the  basic  flow  the  lines  are 


i  -  -1  ,  ft  -  Si 

and  in  the  perturbed  flow 

t  -  -1(1  +  0  R“2+(2/Of) j  >  |  +  6  r‘2+(2/o) j 

We  are  mainly  interested  in  the  value  of  0. 

The  boundary  conditions  at  the  shock  are  obtained  by  solving  the  con¬ 
servation  equations  with  cq  taken  into  account  to  first  order.  With  the 
problem  expressed  in  terms  of  u  and  c,  0  does  not  appear  explicitly 
and  we  have 

u  -  D  U  -  D'  U'1 
-c  -  EU  +  E'd'1 


where 


E'  - 


67-72-l 


2(7+1)  ^27  (7-i) 


D'  - 


2  2 
7+1  Co 


and  U 

-  -  R1"^1/a)  (1  +  0  r“2+<2A*>) 

t 

Hence  at  _  -  -E 

V 

we  have 

r(-E) 

-  -  d'  j 

I  (A. 2) 

s(-E) 

■  ■  ?2a  [■*“  m  +  a  ]  -  «'  ^ 

> 
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where  o  denotes  the  value  at  n  ■  s 


and  the  third  implies 


which  is 


ja(r  s  +r  .  )  +  ( 1-a) (k-1) (r  s  -r  s  ) 

DO  O  O  OOOO 


_  _  (r^-k)b^ 

(r  -s  )  - — — -rrr-  ~  (r  -s  )(r  -ks  )  -  0 

oo  2oo(rj4-k)  o  o  '  o  o 


Hence  T  ,  s  are  known  in  terms  of  6. 
o  o 

Let  us  take  any  trial  value  of  6,  say,  and  suppose  that  the 

values  of  r,  s  at  s  ■  -E  thus  obtained  are  denoted  by  r„^°\  au^°\ 

n  H 

respectively.  Since  the  system  is  linear  any  multiple  of  this  solution  is 
also  a  solution  of  the  differential  equations.  We  require  that  a  multiple 
X  of  this  solution  satisfies  (4.3),  i.e., 


XF  <o) 

XrH 


Ax  P  -  D* 


«„(o>  -  *2 


in  which  X,  P  are  unknown.  Thus  we  can  evaluate  |3  and  also  X  if  it 
is  desired  to  tabulate  the  solution  (this  is  not  done  here). 

The  application  of  Whitham's  rule,  that  the  flow  variables  at  the 


shock,  5  ■  -1,  satisfy  the  characteristic  condition  at  i  ■  gives 

a  simple  algebraic  formula  for  the  approximate  evaluation  of  ^ 

„  D’E-K’P  (l-kHD'E+PE1) 

Pw  2PE  2(D+E) (D+kE) 


However,  one  cannot  expect  this  approximation  to  have  much  accuracy. 
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The  evaluation  of  P  by  this  method  is  extremely  simple  due  to  the 
fact  that  it  is  ontirelv  altrebrnic.  Th»  srWp  nnnl  (oH  tn  eysluats  p 
exactly  (Welsh,  1966)  could  also  be  applied  here.  However,  this  involves 
the  exact  integration  of  the  system  (4.1).  Instead  of  doing  so  it  is  pro¬ 
posed  to  apply  an  approximate  algebraic  method,  the  accuracy  of  which  may 
be  gauged  by  comparing  the  results  of  this  method  to  the  gaB  case  with  the 
exact  values.  In  theory,  one  could  develop  the  exact  regular  solution  of 
(4.1)  as  a  power  series  about  s  -  sq,  although  the  algebra  involved  be¬ 
comes  impossible  after  the  second  terms,  whose  values  are  needed  for  the 
exact  method.  We  shall  make  use  of  the  fact  that  the  range  of  Integration 
is  quite  small,  of  order  0.2,  and  evaluate  the  series  solution  only  as  far 
as  the  second  term,  of  order  s  -  sq.  In  other  words,  we  shall  evaluate 
r^,  sl  for  s  -  and  replace  5^°'  on  the  left  hand  side 

of  (4.5)  by 


'„<0)  -  r„(0)  * 


•hM  ■  +  W*.’ 


and  denote  the  value  of  P  so  obtained  by  P^  (the  first  order  solution). 

The  less  accurate,  zero  order  solution  Pq  will  also  be  found  by  truncating 

the  series  after  a  single  term  so  that  7  ^  are  replaced  by  r 

n  n  O 

s,Q^°\  respectively. 

The  values  of  pQ,  P^,  P^  are  given  for  7  ■  7,  j  -  1,2  In  Table  II 
along  with  the  entropy  varying  case  for  7  *»  3  (Welsh,  1966). 


5.  THE  APPROXIMATE  SOI.UTIOW  TO  THE  COMPLETE  MOTION 

The  work  of  this  section  arises  from  the  work  of  Chester  (1954), 
Chisnell  (1957),  and  Whitham  (1958),  on  plane  shock  waves  propagating  along 
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a  tube  of  varying  cross-sectional  area  A.  Evidently  A  *  R^,  j  *  1,2, 
CCfi'C-pCudS  vG  Wtu^e  «uu  cone -shaped  tubes  and  hence  to  complete  cylindrical 
and  spherical  shock  fronts.  The  method  employed  by  these  authors,  which  is 
in  effect  to  apply  Whitham's  rule  taking  complete  account  of  the* state  ahead 
of  the  shock,  and  so  is  not  restricted  to  strong  shocks,  gives  good  agreement 
with  the  numerical  solution  of  the  complete  system  of  partial  differential 
equations  obtained  by  Payne  (1957). 

The  characteristic  condition  is 


du  -  kdc 


y  cdqj  + 


dR 

R 


where  the  cp  term  is  omitted  for  water, 


2c  , 

7+1  “  M  * 


C 

o 

Cm)M 


n/~C27M2-  (7-1  >  }  ( (7  - 1  )M2+2  ) 


,  since  M  -  jj-  <  0 


This  results  in  the  formula 

JdR  .  2M  dM  - 

ir  +  7^  m  ■  0 

for  the  shock  propagation.  If  we  denote  K(M)  for  a  gas  by  K^°^(M), 
then  Chester's  formula  is 


___2 _ 

K(o)(M) 


2. CM2.-!) 
'/"xy~ 


where  x  -  (7-l)M2+2 
y  -  27M2-7+l 


(5.1) 


and  the  formula  for  water  in  which  9  Is  neglected  is 
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_ 2__ 

K(1>(M) 


2_  7MT+1  I 
M2  J 


(5.2) 


In  the  limit  M-*  1,  for  weak  shocks,  K -*  1/2  in  both  cases  for  all 
7,  so  that  the  integral  for  M  is 

M  -  1  «  R_J/2  , 

the  law  for  acoustic  propagation. 

In  the  limit  ns  M-*  «  the  integral 

K  .  R-U/2K 

so  that  K  "  (l/o  )-l,  spherical  case. 

oo  V 

Graphs  of  K^°^(M)  and  K^(M)  are  given  in  Fig.  2  for  several  cases, 
the  one  most  relevant  to  the  present  work  being  K^(M)  for  7  «  7.  These 
results  give  an  idea  of  the  effect  of  neglecting  entvopy  variation,  which 
is  substantial  for  7  m  1.4,  but  decreases  as  7  increases,  and  is  not 
nearly  so  great  for  7  m  7.  For  a  given  7,  K^(M)  <  K^°^(M)  for  all 

M,  so  that  the  entropy  being  neglected  causes  the  shoe:  to  accelerate  more 
rapidly.  K  is  mainly  monotonic,  increas&&  for  large  7  and  decreasing 
for  small  7.  For  7  m  2,  K  ^  «  1/2  and  K^(M)  varies  very  .. 

so  that 


M 


R'J/2' 


is  a  very  accurate  approximation  for  all  M. 

The  results  also  relate  to  the  propagation  of  a  shock  wave  in  a  tube 
of  varying  cross-section,  the  general  area  charge  dA/A  replacing  the 
term  jdR/R. 
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6.  THE  SPHERICAL  AND  CYLIMDRTCAT.  SHOCK  TL'SZ 

The  self-similar  solution  for  the  final  stages  of  the  collapse  of  «u 
imploding  shock  wave  is  independent  of  the  means  by  which  the  shock  is 
initiated,  except  for  scaling  factors.  This  also  applies  to  perturbations 
to  the  self "Similar  solution  obtained  by  taking  account  of  the  pressure 
ahead  of  the  shock. 

Here  we  shall  consider  the  initial  propagation  of  the  shock  wave  with 
the  initial  state  of  a  spherical  pressure  discontinuity  separating  two  uni¬ 
form  regions.  The  appearance  of  the  initial  radius  of  the  sphere  r^  as 
a  parameter  means  that  the  flow  is  non-similar.  However,  in  the  limit  aa 
the  instant  of  release  is  approached,  the  flow  is  independent  of  and 

is  given  by  the  classical  solution  for  the  one -dimensional  shock  tube, 
caused  by  the  release  of  the  plane  diaphragm  separating  two  uniform  media 
at  different  pressures.  This  solution  is  self-similar.  The  only  parameters 
have  the  dimensions  of  pressure  and  density  (rQ  neglected)  so  that  the 
similarity  variable  is  r/t  and  the  flow  is  uniformly  expanding.  In  this 
solution  the  pressure  ahead  of  the  shock  is  taken  into  account  exactly 
(for  the  final  collapse  the  similarity  variable  is  r/t  ,  0  <  a  <  1  and 
similarity  demands  that  the  pressure  ahead  of  the  shock  be  neglected). 

As  the  shock  propagates  initially  dlaturbances  will  arise  in  the  flow 
due  to  the  departure  from  plane  geometry,  and  these  will  cause  an  imploding 
shock  to  accelerate  and  an  expanding  one  to  decelerate.  These  disturbances 
will  be  treated  as  a  perturbation  on  the  basic,  one -dimensional,  self¬ 
similar  solution,  to  obtain  the  initial  acceleration  of  the  shock  wave. 

The  analysis  will  be  applicable  to  any  combination  of  gas  and  water  for 
the  two  media,  and  al6o  to  either  the  sphere  being  at  a  lower  pressure  than 
its  exterior  or  the  reverse  situation.  Since  it  is  known  that  the  final 


collapse  rate  is  singular,  a  aorioa  obtained  a*  a  n»rtiirh«H™  “ill  cer¬ 
tainly  not  be  valid  up  to  the  moment  of  collapse,  and  can  only  describe 

the  initial  motion  (i-e.,  for  small  values  of  [r-r  J/r  ). 

o  o 

The  representation  of  the  one-dimensional  shock  tube  flow,  with  r  -  tQ 
as  the  diaphragm  position  is  shown  in  Fig.  4.  For  this  diagram  the  high 
pressure  region  is  r  <  rQ.  The  gas  in  the  lower  pressure  region  is  com¬ 
pressed  through  a  shock  wave  and  that  in  the  high  pressure  region  is 
expanded  through  a  dimple  wave.  The  regions  0,4  beyond  the  influence  of 
these  waves  are  in  the  initial  stationary  states.  The  regions  2,3  are 
in  uniform  motion  with  the  same  pressure  and  fluid  velocity.  The  contact 
front  separates  the  two  media.  The  region  1  is  a  point -centered  simple 
rarefaction  wave. 

In  the  one-dimensional  system  this  flow  would  persist  for  all  time  in 
the  absence  of  any  dissipation.  The  disturbances  on  this  flow  due  to  a 
spherical  or  cylindrical  geometry  will  start  from  t  ■  0  and  will  be  small 
in  the  initial  stages.  These  disturbances  will  travel  along  characteristics, 
so  let  us  consider  the  overall  picture  of  the  disturbances.  Consider  first 
the  fluid  in  the  high  pressure  region  r  <  rQ,  i.e.,  the  regions  0,  1,  2. 

No  disturbances  can  propagate  in  the  region  0,  which  remains  uniform. 

Also  the  path  of  the  last  characteristic  of  the  simple  wave,  separating  0 
and  1,  is  not  disturbed.  In  the  expansion  wave,  region  1,  there  will  be 
one  degree  of  freedom  for  each  family  of  characteristics  along  which  dis¬ 
turbances  may  propagate.  By  a  degree  of  freedom  we  mean  an  arbitrary 
constant  in  the  general  solution  for  the  region,  and  these  constants  which 
arise  will  be  determined  from  appropriate  boundary  conditions.  In  any 
region  there  are  three  families  of  characteristics,  c^  and  c^ 
defined  by  d/dt  *010,  d/dt  «  u,  respectively.  Thus  in  any  region 
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there  will  be  at  most  three  degrees  of  freedom.  There  can  be  no  entropy 
disturbance  in  the  simple  wave,  so  that  there  is  no  degree  of  freedom 
associated  with  the  c''"  characteristics,  i.e.,  the  particle  paths. 

The  negative  characteristics  are  intersecting  straight  lines  through 
r  «  r  ,  t  *  0.  In  the  basic  flow  this  point  is  singular  as  the  character¬ 
istics  intersecting  here  have  different  values  of  the  characteristic  vari¬ 
able.  Since  all  perturbations  must  be  zero  at  t  *  0  and  all  these  lines 
intersect  at  t  ■  0,  it  follows  that  there  can  be  no  disturbances  propa¬ 
gating  along  Hence  there  is  only  one  degree  of  freedom  in  the 

f+) 

simple  wave  due  to  disturbances  along  c  ,  and  this  will  be  determined 
by  the  condition  that  it  is  zero  on  the  last  characteristic  of  the  simple 
wave. 

There  is  no  entropy  disturbance  in  region  2  but  there  will  be  disturb¬ 
ances  on  c<+)  and  c^  \  the  former  coming  from  the  simple  wave  and  the 
latter  from  the  shock.  The  solution  in  this  region  has  to  be  matched  with 
that  on  the  first  characteristic  of  the  simple  wave.  This  condition  deter¬ 
mines  the  constant  associated  with  c^  in  2  to  that  of  c^  In  1 
as  the  disturbances  propagating  along  c^  ^  in  region  2  do  not  enter 
into  this  matching. 

Thus  the  motion  of  the  gas  to  the  left  of  the  contact  front  can  be 
described  except  for  the  degree  of  freedom  associated  with  ^  in  region 
2,  which  will  be  determined  by  matching  this  solution  with  that  for  the 
gas  initially  in  r  >  rg,  along  the  contact  frdnt. 

The  region  4  ahead  of  the  shock  is  unperturbed.  However,  the  path 
of  the  shock  itself  will  be  altered,  causing  disturbances  along  all  three 
families  of  characteristics  in  the  region  3  behind  it.  The  conservation 
equations  across  the  shock  serve  to  determine  the  boundary  values  of  the 
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solution  in  terras  of  the  change  in  shock  speed,  so  that  the  flow  in  this 
region  is  expressible  in  terms  ot  one  unknown,  the  change  in  shock  speed. 

There  arc  two  conditions  to  be  satisfied  at  the  contact  front,  the 
continuity  of  pressure  and  fluid  velocity.  Since  the  boundary  valueB  to 
the  left  and  right  of  the  front  each  contain  one  unknown,  these  two  con- 
ditions  will  determine  the  flow  completely. 

The  characteristics  along  which  disturbances  propagate  are  sketched 
in  Fig.  3.  The  problem  is  greatly  simplified  if  the  disturbances  propagating 
along  c  in  region  3  are  neglected,  as  the  propagation  of  the  shock 
wave  is  then  determined  simply  by  examination  of  the  boundary  values  at  the 
shock,  without  reference  to  the  motion  behind  the  shock.  This  situation 
would  arise  if  there  were  no  simple  wave  and  the  regions  0,1  were  identi¬ 
cal  to  2.  An  extensive  region  of  uniform  flow  behind  the  Bhock  cannot  arise 
in  an  entirely  spherical  or  cylindrical  geometry.  However,  it  describes  a 
uniform  plane  shock  wave  traveling  in  a  one -dimensional  tube  and  meeting  an 
appropriately  varying  area  change.  This  is  the  case  studied  by  Chester 
(19^4),  Chlsnell  (1957),  and  Gundersen  (1958).  For  a  very  weak  shock,  which 

(4.) 

is  nearly  sonic,  the  disturbances  reaching  it  along  c  take  a  long  time 
to  do  so  and  hence  their  effect  will  be  increasingly  small  for  increasingly 
weak  shocks.  The  simplest  derivation  of  the  formula  for  the  shock  propaga¬ 
tion  neglecting  these  disturbances  is  that  of  Whitham  (1958),  and  it  will 
be  shown  that  neglecting  them  in  the  present  analysis  reduces  to  his  formula. 
The  release  of  a  pressurized  sphere  in  air  has  been  studied  by  Friedman  (1961), 
who  uses  Whitham's  rule  for  the  propagation  of  the  main  shock  wave. 

Equationo  of  Motion 

The  equations  of  motion  for  the  spherically  or  cylindrically  symmetric 
flow,  in  the  absence  of  dissipation  due  to  heac  conduction,  are  those  given 
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previously  in  Section  3. 

a*  * 

these  are 

«t*  +  (uic)^* 

<Pt  +  “Pr  -  0 

The  basic  self-similar  flow  Is  given  by  the  above  equations  with  the  geometry 
term  Involving  j  omitted,  the  variables  being  functions  of  (r-rQ)/t  only. 
The  perturbations  to  be  evaluated  here  are  essentially  due  to  the  effect  of 
this  term. 

The  basic  solution  is  of  the  form 


Defining  the  characteristic  variables 


—  ti  +  — i—  /. 

2  "  “  r-1  ” 


-  j  12 

-  +  27  uc  +  27  C  Vx 


(6.1) 


o>> 


*0  “  <P0<*->  .  X  - 

and  also  UQ  -  constant.  Let  us  expand  about  this  solution  in  the  form 


0*  -  aoa(\)  +  ectjV)  +  e2a2+(\)  +  .  .  .  . 
r-r 

where  e  ■  - £ 

r 


(6.2) 


We  shall  be  concerned  with  only  the  first  order  corrections,  i.e.,  those  with 
suffix  1.  On  substituting  the  variable,  expanded  in  this  form,  into  the 
equations  of  motion  and  equating  like  powers  of  X,  we  obtain  the  equations 
governing  the  basic  flow 


(X-u  )m  '  m 
0^0 

(u  ±c  -\)a  * 

O  0  O 


(6.3) 


*• 

. -/ 


0 
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and  the  equations  for  the  perturbations 

\(u  -X,)«p..  *  +  u  ?,  +  Xu.?,.'  ■  0 

O  I.  Cl  a  w 

*<wMai1'  +x>“i±ci)v' +  iv'/i’ 

■  ‘iV.'b  |2XciV  + 

The  simplicity  of  the  equations  (6.3)  is  due  to  the  combination  of 

self-similarity  and  characteristic  propagation.  The  first  ia  the  statement 

<p  is  constant  except  for  a  possible  discontinuity  on  X  ■  u  ,  i.e.,  the 

o  o 

particle  path  through  the  point  of  initiation  in  the  r-t  plane.  Thus  the 
right  hand  side  of  the  second  ia  zero,  and  this  becomes 

+ » 

(u  ±c  -X)CK  ■  0 

0  0  o 

which  atotea  that  0^*  is  constant  except  possibly  on  an  appropriate  char¬ 
acteristic  through  the  point  of  initiation.  Here,  in  the  case  of  the  point- 
centered  simple  wave,  it  is  possible  for  there  to  be  an  intersecting  family 
of  such  characteristics,  giving  a  region  in  which  the  first  factor  equated 
to  zero  is  the  appropriate  solution  of  this  *»quation.  For  all  other  regions, 
which  are  uniform,  o.*  ■  constant  is  the  solution.  The  characteristics  of 
the  system  (6.4)  are  evidently  the  three  characteristics  of  the  flow. 

The  Simple  Wave 

Here  we  shall  solve  the  system  (6.3)  for  the  familiar  solution  for  a 
point-centered  simple  wave  and  hence  the  system  (6.4)  for  the  perturbations. 
In  these  equations  a  double  sign  ±  was  used  for  simplicity,  the  upper 
sign  referring  to  the  c^  equation  and  the  lower  to  the  c^’\  In  order 
to  develop  the  two  cases  of  high  pressure  inside  and  outside  the  sphere 
simultaneously  a  similar  notation  will  be  used.  The  former  will  be  referred 
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to  as  case  +,  and  will  be  given  by  the  upper  signs  in  the  eauationa. 

the  later,  case  will  be  given  by  the  lower  signs.  Thus  the  two  signs 

refer  here  to  different  physical  situations,  whereas,  as  used  previously, 

they  refer  to  different  equations  describing  the  same  physical  situation. 

The  basic  simple  wave  is 


u  +  c  -  X 
o  o 


The  solution  for  is 

*1  -  r  ucor  +  O'- dm' 

and  the  equations  for  are 


7+1 
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In  each  pair  of  equations  one  is  algebraic  so  that  their  solution  has  only 
one  arbitrary  constant.  This  is  due  to  there  being  no  disturgances  propa¬ 
gating  along  the  rectilinear.  Intersecting  characteristics,  the  similarity 
approach  ruling  out  such  disturbances  automatically. 

The  solution  of  the  differential  equations  which  appear  in  (6.6)  and 
(6.7)  may  be  expressed  as 

_ -J  -  J  c* 


(7-1X37-5) 

*  _ 6 _ 

27(37-1) 


(+}  - 


'oR 
(7-1) (3-7) 


{+) 


(+) 


2Z_ 

7-1 


+  Kj  (+) 


7+1 

2(7-1) 


(6.8) 


where  {+}  -  (2cqR  +  (7-l)X) 


7/3,  5/3  , 


and  K.  are  arbitrary  constants. 


The  algebraic  equations  may  be  combined  as 

x  [  "l  +  r^r  ClJ  '  1  A  <^lcoR>  <+> 


(6.9) 


Since  the  perturbed  flaw  is  homentropic  A  may  be  set  equal  to  zero 
in  (6.8).  The  flow  is  continuous  across  the  last  characteristic  of  the 
simple  wave  separating  the  regions  0  and  1,  i.e.,  X.  -  +  coR,  so  that 

u^,  c^  must  be  zero  on  this  line.  The  equations  (6.8),  (6.9)  for  u^,  c^ 
are  singular  only  In  the  case  7  »  1.  The  right  hand  side  of  (6.9)  is  zero 
on  \  ■  +  ceR,  so  K*  is  determined  by  equating  the  right  hand  side  of 
(6.8)  to  zero.  This  gives 


37-5 


2  j  (7+1) 


2(7-1)  „  2(7-1) 

coR 


(37-5)(3-7) 


(6.10) 


The  solution  in  the  simple  wave  is  now  determined  completely  by  (6.8),  (6.9), 


with  A  «  0  and  K  given  by  (6.10). 
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In  the  basic  simple  wtve  the  first  characteristic  is  the  line  on  wnicr. 

the  value  of  the  fluid  ve  loci  tv  is  the  same  as  that  behind  the  shock*  u  • 

r 

say.  On  this  line 


u 

o 


u 

8 


"  coR  +  2TIub 

7+1  - 

m  *  -  —  U  4*  C  D 
2  6  oR 


The  Uniform  Regions  2  and  3 

In  the  case  of  a  uniform  homent.ropic  basic  flow  the  equations  (6.4) 
governing  the  perturbations  are 


\(uQ-\>V  +  uo*l 


\(uo±co-\)aia'  +  <uo±co)a1±  -  +  f  Vo  +  #"  ^l'**!* 


where  u  ,  c  are  constant, 
o’  o 

The  solution  of  this  system  is 


\-u 


o 

.  Ac 

^i*  ■  ^iuoco  +  ir+ K'±(uo±coA) 


(6.11) 


The  last  term  is  constant  on  uo±cq-\  «  X,  i.e.,  on  c^,  but  the  other 
t^rms  are  not  associated,  with  characteristic  propagation  as  written.  The 


entropy  term  is  associated  with  c 
2 


(o) 


■o  we  can  write  it  as 


Ac 
_ o_ 
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-  *  z?  <N-u0)  ±  jT  <Uo±Co“X> 


and  the  first  term  may  be  written 

) 


,  u  c 
2  o  o 


»  u  c  \ 
J  op 


.  u  ±c 

,  J  a  o 

+  •=  u  c  - 


2  u  ±c  2  o  o  u  ±c 
O  O  0  0 


•  A- 


so  that  (6.11)  may  bo  expressed  In  the  form 
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Xa, 


-  +  rr 


\  u  c  X  Ac  ,  \ 

'  ±  sr2  (X-u  )  +  K^(u  ±c  -X)  / 

o  '  o  o  / 


2  u  ±e 
o  o 


rr 


(6.12) 


where  in  the  last  equation  the  terms  are  respectively  the  geometry  modifi¬ 
cation,  the  entropy  disturbance  along  c^°\  and  the  disturbance  traveling 
along  c(±*.  By  neglecting  the  disturbances  reaching  the  shock  from  behind, 
l.e.,  by  setting  the  appropriate  K  zero,  Whifchara's  formula  is  obtained, 

A  retained  for  a  gas  and  omitted  for  water. 

The  particles  in  the  region  2  between  the  simple  wave  and  the  contact 
front  do  not  pass  through  the  shock,  so  that  A  ■  0  in  region  2,  and 

*“1*  -  -j#+ V<VV*> 

o  o 

where  K^±  are  both  non-zevo. 

This  flow  has  to  match  identically  with  the  perturbed  flow  in  the  simple 
wave  on  the  first  characteristic  of  the  simple  wave  separating  the  two 
regions,  which  determines  one  of  .  The  remaining  Kj  will  be  deter¬ 
mined  by  the  conditions  at  the  contact  front.  Here 


so  that 


,  1  u  c  c  , 

a*  -  +  J-  -2~-  ±  -a  K  * 

1  2  u  ±c  u  2 

o  o  o 


(6.13) 


where  u 


X  2x1 

CoR  +  2  ut 


The  solution  in  the  region  3  between  the  contact  front  and  the  shock 


is 
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Xtpj  -  A(X.-uq) 


to.!*; 


i  u  c  k  Ac  . 

+  4-  -2-2-  ±  tt-2-  (A- u  >  +  K.„  (u  ±c  -k) 

i  u  ±c  tr  O'  j  '  o  o 

o  o 


If  the  medium  in  this  region  is  gas  then  A  is  non-zero  and  is  determined 
by  the  conservation  equations  at  the  shock.  For  water,  in  which  the  effect 
of  entropy  variations  is  neglected,  A  will  be  zero.  The  conservation 
equations  at  the  shock,  i.e.,  on  X  -  Uq,  give  three  equations  relating 
A,  K^1,  and  U1#  so  that  the  solution  in  this  region  is  expressible  in 
terms  of  the  single  unknown  U^.  The  boundary  values  at  the  contact  front 
are  given  by  (6.14)  with  k  ■  uq.  The  physical  situations  are  here  dis¬ 
tinguished  by  the  fact  the  uq,  Uq  are  negative  for  the  imploding  shock 
and  positive  for  the  exploding  shock  (cq  being  always  positive). 

The  continuity  of  fluid  velocity  and  pressure  at  the  contact  front 
give  two  equations  involving  the  unknown  and  U^.  No  calculations 

have  been  performed  as  yet. 
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